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Motivated by the ongoing debate about nanophotonic control of Forster resonance energy transfer 
(FRET), notably by the local density of optical states (LDOS), we study an analytic model system 
wherein a pair of ideal dipole emitters - donor and acceptor - exhibit energy transfer in the vicinity of 
an ideal mirror. The FRET rate is controlled by the mirror up to distances comparable to the donor- 
acceptor distance, that is, the few-nanometer range. For vanishing distance, we find a complete 
inhibition or a four-fold enhancement, depending on dipole orientation. For mirror distances on 
the wavelength scale, where the well-known ‘Drexhage’ modification of the spontaneous-emission 
rate occurs, the FRET rate is constant. Hence there is no correlation between the Forster (or 
total) energy transfer rate and the LDOS. At any distance to the mirror, the total energy transfer 
between a closely-spaced donor and acceptor is dominated by Forster transfer, i.e., by the static 
dipole-dipole interaction that yields the characteristic inverse-sixth-power donor-acceptor distance in 
homogeneous media. Generalizing to arbitrary inhomogeneous media with weak dispersion and weak 
absorption in the frequency overlap range of donor and acceptor, we derive two main theoretical 
results. Firstly, the spatial dependence of the Forster energy transfer rate does not depend on 
frequency, hence not on the LDOS. Secondly the FRET rate is expressed as a frequency integral of the 
imaginary part of the Green function. This leads to an approximate FRET rate in terms of the LDOS 
integrated over a huge bandwidth from zero frequency to about 10 x the donor emission frequency, 
corresponding to the vacuum-ultraviolet. Even then, the broadband LDOS hardly contributes to the 
energy transfer rates. Using our analytical expressions, we plot transfer rates at an experimentally 
relevant emission wavelength A = 628 nm that reveal nm-ranged distances, and discuss practical 
consequences including quantum information processing. 
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ratus of plants and bacteria 00. Many applications 
are based on FRET, ranging from photovoltaics 0, [Hi, 
lighting 0,1, and magneto-optics :9], to sensing 0 
where molecular distances [TT|, [l^] , and interactions are 
probed 00 . FRET is also relevant to the manipula¬ 
tion, storage, and transfer of quantum information 0 - 
0 - 

Modern nanofabrication techniques have stimulated 
the relevant question whether Forster transfer can be 
controlled purely by means of the nanophotonic environ¬ 
ment, while leaving the FRET pair geometrically and 
chemically unchanged. In many situations, the effect of 
the nanophotonic environment can be expressed in terms 
of the local density of optical states (LDOS) that counts 
the number of photon modes available for emission, and is 
interpreted as the density of vacuum fluctuations [21 1. 122 . 
An assumption behind many recent FRET studies has 
therefore been that if there is an effect of the nanopho¬ 
tonic environment on FRET rates, then it should be pos¬ 
sible to find a general law that describes the functional 

_ dependence of FRET rates on the LDOS. While it is an 

assumption, it is a fruitful one as it allows for experi- 
* Electronic address: mwubs@fotonik.dtu.dk mental verification. Curiously, different dependencies of 

t Electronic address: w.l.vos@utwente.nl, world wide web: www. photon it IKifcclfta |bnl:OHi On the LDOS were reported in a number 


I. INTRODUCTION 


A well-known optical interaction between pairs of 
quantum emitters - such as excited atoms, ions, 
molecules, or quantum dots - is Forster resonance en¬ 
ergy transfer (FRET). In this process, first identified in 
a seminal 1948 paper by Forster, one quantum of exci¬ 
tation energy is transferred from a first emitter, called a 
donor, to a second emitter that is referred to as an ac¬ 
ceptor 0. FRET is the dominant energy transfer mecha¬ 
nism between emitters in nanometer proximity, since the 
rate has a characteristic (rp/r^a.) 6 distance dependence 
(with rp the Forster radius and ra a the distance between 
donor and acceptor). Other means to control a FRET 
system are traditionally the spectral properties of the 
coupled emitters - the overlap between the donor’s emis¬ 
sion spectrum and the acceptor’s absorptions spectrum - 
or the relative orientations of the dipole moments 00 - 
FRET plays a central role in the photosynthetic appa- 
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of experimental studies, leading to the debate how the 
Forster energy transfer rate depends on the LDOS. Pio¬ 
neering work by Andrew and Barnes indicated that the 
transfer rate depends linearly on the donor decay rate 
and thus the LDOS at the donor emission frequency [23| . 
as was confirmed elsewhere 2M • The linear relation be¬ 
tween the two rates found in [231] was supported by the 
theory of Ref. [2|| but only fortuitously within a limited 
parameter regime. Subsequent experiments suggested a 
transfer rate independent of the LDOS [27j, a dependence 
on the LDOS squared [26|, or qualitative effects [28l.[29[. 
Possible reasons for the disparity in these observations in¬ 
clude insufficient control on the donor-acceptor distance, 
on incomplete pairing of every donor to only one accep¬ 
tor, or on cross-talk between neighboring donor-acceptor 
pairs. 

Recently, the relation between Forster transfer and the 
LDOS was studied using precisely-defined, isolated, and 
efficient donor-acceptor pairs [3Cj. The distance between 
donor and acceptor molecules was fixed by covalently 
binding them to the opposite ends of a 15 basepair long 
double-stranded DNA. A precise control over the LDOS 
was realized by positioning the donor-acceptor pairs at 
well-defined distances to a metallic mirror [22], [3l}, [32]]. 
The outcome of this experimental study was that the 
Forster transfer rate is independent on the optical LDOS, 
as was confirmed by theoretical considerations. Con¬ 
sequently, the Forster transfer efficiency is greatest for 
a vanishing LDOS, hence in a 3D photonic band gap 
crystal [Hj]. Similar results were obtained with different 
light sources (rare-earth ions), and with different cavi¬ 
ties [HI, [35]. On the other hand, a linear relation between 
LDOS and FRET rate was reported in experiments with 
donors and acceptors at a few nanometers from metal 
surfaces 0, HtJ , while Ref. [11] reported no general the¬ 
oretical relationship between LDOS and FRET rate near 
a metallic sphere. In Ref. [H the measured dependence 
of the FRET rate on the LDOS was reported to be weak 
for single FRET pairs; an observed drop of the total en¬ 
ergy transfer rate of a donor close to a surface was mainly 
attributed to the simple fact that fewer statistically dis¬ 
tributed acceptors are available close to the surface; re¬ 
cent theoretical work on collective e nerg y transfer sup¬ 
ports these results in the dilute limit [40j] . 

Most theoretical papers agree that both the energy 
transfer rate and the spontaneous-emission rate can be 
expressed in terms of the Green function of the nanopho¬ 
tonic medium. One may argue that energy transfer and 
optical LDOS are therefore related. But one may also 
argue to the contrary, since the energy transfer rate de¬ 
pends on the total Green function describing propaga¬ 
tion from donor to acceptor, whereas the spontaneous- 
emission rate depends on the imaginary part of the Green 
function at the donor position only |25j. It is not clear 
whether this situation entails correlations between the 
two quantities and if so, what is their functional rela¬ 
tionship. Therefore, and in view of the different results 
in the literature, we decided that a study of a simple 
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Figure 1: (Color online) We study pairs of donor and acceptor 
dipoles that are separated by a distance rd a , and located at a 
distance 2 : from an ideal mirror. We focus on two configura¬ 
tions, both with dipoles oriented perpendicular to the position 
difference vector of donor and acceptor (/td,/t a ) -J- (ra — r a ): 
(a) Both dipole moments of donor and acceptor are parallel 
to the mirror surface (‘parallel configuration’, ||) and parallel 
to each other; (b) Both dipole moments of donor and accep¬ 
tor are perpendicular to the mirror surface (‘perpendicular 
configuration’, _L) and parallel to each other. 

analytical model is timely. 

In this paper, we first study energy transfer in a pro¬ 
totypical nanophotonic medium, namely the Drexhage 
geometry [3l| . near an ideal mirror. This is one of 
the simplest inhomogeneous dielectric media for which 
position-dependent spontaneous-emission rates are ana¬ 
lytically known [32]. Here we show that the spatial de¬ 
pendence of both the total and the Forster resonant en¬ 
ergy transfer rates can be calculated analytically. Such 
exact results have a value of their own, and allow for a 
critical and straightforward assessment of possible cor¬ 
relations between the FRET rate and the LDOS. After 
studying the phenomenology near an ideal mirror, we de¬ 
rive general results for arbitrary inhomogeneous weakly 
dispersive media, and thereby find Forster transfer rates 
that generalize the well-known l/(n 4 r® a ) dependence of 
the homogeneous nondispersive medium with refractive 
index n. Specifically, using the energy-transfer theory 
by Dung, Knoll, and Welsc h_[25| and the Green-function 
properties derived in Ref. [4l| as starting points, we de¬ 
rive a simple and important new consequence: for the 
large class of photonic media that have little dispersion 
and absorption in the donor-acceptor frequency overlap 
range, there is no correlation between the LDOS and the 
position-dependent FRET rate. Nevertheless, the FRET 
rate is controlled by the distance to the mirror, but only 
at mirror distances comparable to or smaller than the 
donor-acceptor separation. 

II. PHYSICAL PROCESSES AND GEOMETRY 

We study energy transfer from a single donor to a sin¬ 
gle acceptor separated by a distance ru a = | r a — r t j |. To 
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limit parameter space, we focus on situations in which 
the donor and the acceptor have the same distance z to 
the mirror, and where the dipole moments of dipole and 
acceptor point in the same direction. In the parallel (||) 
configuration shown in Figure QJa,), both dipole moments 
are oriented parallel to the mirror, and the dipoles point 
normal to the mirror in the perpendicular (_L) configura¬ 
tion of Figure IHb). 

The total energy transfer rate 7d a between a donor and 
an acceptor dipole in any nanophotonic environment is 
given by 

/ OO 

dwcra(w) w(r a ,rd,w)<7d(u;), (1) 

-OO 

where Cd, a (w) are the donor (single-photon) emission 
and acceptor (single-photon) absorption spectra in free 
space [25|, [42|. All effects of the nanophotonic environ¬ 
ment are contained in the transfer amplitude squared 
w( r a ,rd,w) that can be expressed in terms of the Green 
function G(r a ,rd,w) of the medium, and the donor and 
acceptor dipole moments fi diMa respectively, as 

w(r a ,r d ,w) = ^ j |/4 ' G(r a ,r d ,a;) • /r d | 2 - (2) 

These expressions for the total energy transfer rate were 
originally derived in an important paper by Dung, Knoll, 
and Welsch for a general class of nanophotonic media 
that may exhibit both frequency-dispersion and absorp¬ 
tion [25j [57l ]. The total energy transfer rate is the com¬ 
bined effect of both radiative and Forster energy transfer 
processes. We will first study total energy transfer rates, 
and in Sec. m we discuss what fraction of this energy 
transfer is Forster resonance energy transfer. 

For the energy transfer rate Eq. © we only need to 
know the Green function in the frequency interval where 
the donor and acceptor spectra overlap appreciably. In 
case of molecules that have among the broadest band- 
widths, this overlap has a typical relative bandwidth 
of only a few percent. Hence it is reasonable to ne¬ 
glect absorption and material dispersion in this narrow 
overlap region. Thus s(r,ui) can be approximated by 
a real-valued frequency-independent dielectric function 
e(r). The corresponding Green function G(r,r',w) is the 
solution of the usual wave equation for light 


Glossary of transfer and emission rates 

^7da 

total donor-accept or energy transfer rate, Eq. Q 

'yda 

narrowband approximation of transfer rate, Eq. @ 

'T’se 

spontaneous emission rate of the donor, Eq. © 

7f 

exact FRET rate from donor to acceptor, Eq. (12111 


broadband LDOS approximated FRET rate, Eq. (12711 

7f P) 

high-frequency approximated FRET rate, Eq. (12911 


Table I: Symbols for the various energy transfer and emission 
rates used in this paper, with their defining equations. 


can approximate the energy transfer rate by 

/ OO 

dw(7 a (w)(7 d (w), (4) 

-OO 

where cu ^ a is the frequency where the integrand in the 
overlap integral assumes its maximal value. The overlap 
integral is the same for all nanophotonic environments, 
so that the ratio of energy transfer rates in two different 
environments only depends on the ratio of w( r a , r d ,w da ) 
in both environments. 

Spontaneous emission of the donor is a process that 
competes with the energy transfer to the acceptor. The 
donor spontaneous-emission rate 7 se (r, D) at position r 
with real-valued dipole moment fi = /ifi and transition 
frequency w d is expressed in terms of the imaginary part 
of the Green function of the medium as 

^ H ■ Im[G(r d , r d ,Wd)] • M (5) 

or 7 (r d ,w d ,/4 = 7r/i 2 w d p p (r d , w d , /x)/(3fe 0 ) in terms 
of the partial LDOS p p ( r d ,w d ,/i) = — (6w d /7rc 2 )/i • 
Im[G(ra, r d , a>d)1 • A, where jj, is a dipole-orientation unit 
vector [2l], [43|. The optical density of states (LDOS) 
is then defined as the dipole-orientation-averaged partial 
LDOS [43}. In general both the LDOS and the partial 
LDOS for any dipole orientation are fixed once the partial 
LDOS is known for nine independent dipole orientations, 
but for planar systems considered here, the two directions 
T and || suffice for a complete description [44J [S9| . We 
do not average over dipole orientations, as we are inter¬ 
ested in possible correlations between energy transfer and 
spontaneous-emission rates for a fixed dipole orientation. 


7se(r d ,Wd) = 


-(■ 


2 7i 


\fl£ 0 C 2 


— V x V x G(r, r', oj) +e(r) ^ G(r,r / ,w) = <5(r — r')l, 

_ (3) 

with a localized source on the right-hand side |58|. Un¬ 
like e(r), the Green function G(r, r',w) is frequency- 
dependent and complex-valued. 

While the energy transfer rate in Eq. © evidently 
depends on the donor and acceptor spectra (7 d (w) and 
<T a (w), we are in this paper more interested in the de¬ 
pendence on the environment as given in Eq. ©. We 
therefore assume that the donor and acceptor overlap in 
a narrow-frequency region in which the transfer ampli¬ 
tude w(r a ,r d ,u;) varies negligibly with frequency, so we 


To proceed we need to compute Green functions. First, 
the Green tensor in a homogeneous medium with real¬ 
valued refractive index n is given by [4^| 

G h (ri,r 2 ,w) = G h (r,w) = 

e w 1 

— A — [P(w)\ + Q(w)r <g>r] + —-y-^<5(r)l, ( 6 ) 

47rr 3 (nw/c) 

where r = 17 — r 2 , the functions P,Q are defined as 
P(w) = (1— ■u> -1 +u’ -2 ) and Q(w) = (—1+31W 1 — 3iW 2 ), 
and the argument equals w = (incur/c). For n = 1, 
Gh equals the free-space Green function, denoted by Go. 
For distances much smaller than an optical wavelength 
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(r = |r| <C A = 2nc/(nu})), the Green function scales as 
G h (r,w) oc l/(n 2 r 3 ). From Eqs. © and (0 we then ob¬ 
tain the characteristic scaling of the Forster transfer rate 
as 7da oc l/( n4r da) : the Forster transfer rate strongly 
decreases with increasing donor-acceptor distance and 
with increasing refractive index. In contrast, it follows 
from Eq. © that the spontaneous-emission rate 7 se in 
a homogeneous medium is enhanced by a factor n com¬ 
pared to free space. More refined analyses that include 
local-field effects likewise predict a spontaneous-emission 
enhancement [47| . These major differences between the 
energy transfer rate and the emission rate in a homoge¬ 
neous medium already give an inkling on the behavior in 
nanophotonics. 

Next, we determine the Green function of an ideal flat 
mirror within an otherwise homogeneous medium with 
refractive index n. While the function can be found 
with various methods f d2i . |50], we briefly show how it 
is obtained by generalizing the multiple-scattering for¬ 
malism of Ref. |48| for infinitely thin planes. In the usual 
mixed Fourier-real-space representation (k|| , z) relevant 
to planar systems with translational invariance in the 
(a;, y)-directions, the homogeneous-medium Green func¬ 
tion Gh(kii, z, z', u>) becomes 


G h = 


10 0 

0 k 2 z -k\\k z s zz ' 
0 -fc|| k z s zz > /cjj 



5{z-z’) „ 

(' nuj/c ) 2 


(7) 


where the scalar Green function is given by g^ = 
5 h(k||, 2 , 2 ',ci;) = exp(2ik z \z-z'\)/(2ik z ), k z = ( nco 2 /c 2 - 
fc|) 1//2 , s zz i = sign(z — z') and the matrix is represented 
in the basis (sk,Pk,z), where k is the wave vector of 
the incoming light, z is the positive-z-direction, §k is the 
direction of s-polarized light (out of the plane of inci¬ 
dence), and pk points perpendicular to z in the plane of 
incidence. An infinitely thin plane at z = 0 that scatters 
light can be described by a T-matrix T(k||,w), in terms 
of which the Green function becomes 


functions as 


G (z,z') = G^{z-z')-Gh{z + z')+2 ) g h (z+z')zz, 

\nu) J 

(9) 

where the (k| , w)-dependence of the Green functions was 
again suppressed. To understand energy transfer rates 
near a mirror, we need to determine the Green function 
in the real-space representation, which is related to the 
previous equation by the inverse Fourier transform 

G(r,r',u;) = J d 2 k||G(k||,z,z',u;)e* k ii’ (p_p,) , 

(10) 

where p = (x, y) and p' = (x',y') so that r = ( p,z ). 
Knowing that this inverse Fourier transform when ap¬ 
plied to Gh(k||,z — z',uj) leads to the expression © 
also helps to evaluate the transform of Gh(k 11 , z + z', w). 
Similarly, the third term on the right-hand side of 
Eq. © transform analogous to the zz-component of the 
homogeneous-medium Green tensor. We thus find the 
Green function for an ideal mirror within a homogeneous 
medium as the sum of three terms: 

G(r,r',w) = G h (r, r', w) — G h (p, 2 + z' 1 p', 0, u>) 

+2Gq Z (p, z + z', p', 0,w)zz. (11) 

For the parallel configuration, we find 

e inur da / c 

P ■ G(r a , r d , ca) • /x = -p — - P{inujr da /c) 

47rr da 

^inuju/c 

+ /x 2 — - P(inum/c ), (12) 

47 TU 

where r da is the donor-acceptor distance, z the distance 
of both dipoles to the mirror, and u = [r da + ( 22 :) 2 ] 1 / 2 , 
and /xll = py as in Fig. [1] Inserting this result into 
Eq. © immediately gives the squared transfer amplitude 
u>(r a ,r d ,w) of Eq. © for the parallel configuration. 

For the perpendicular configuration we find 


G (z, z') = G h (2, z) + G h (z, 0)T G h (0, z'), (8) 

where the (k||, cu) dependence was dropped. It was found 
in Ref. 0 that for an infinitely thin plane that models 
a finite-thickness dielectric slab of dielectric constant e, 
the T-matrix assumes a diagonal form in the same basis 
as Gh in Eq. ©, in particular T = diag(T ss , T pp , 0). The 
infinitely thin plane becomes a perfectly reflecting mir¬ 
ror if we choose for example a lossless Drude response 
with e = l—io 2 /ijj 2 , in the limit of an infinite plasma fre¬ 
quency uj p —> oo. Hence the T-matrix for a perfect mirror 
in a homogeneous dielectric has nonzero diagonal compo¬ 
nents T ss = —2 ik z and T pp = —2 i(nui/c) 2 /k z . The ideal 
mirror divides space into two optically disconnected half 
spaces, and below we only consider the half space z > 0. 
It then follows that the Green function for the ideal mir¬ 
ror is written in terms of homogeneous-medium Green 


ginojrda/c 

/x k • G(r a ,r d ,w) • p 1 - = -p 2 — -P(inwr da /(|3) 


/<= 


n inuju/c 


47TU 


47T r da 
2 


P{inuju/c) + 4 ^ Q(inoju/c ) 


with /x k = pz as in Fig. |T|and u as i n Eq- m, whereby 
the squared transfer amplitude of Eq. © is also deter¬ 
mined for the perpendicular configuration. 

For completeness, we also give single-emitter 
spontaneous-emission rates near the mirror (neglecting 
local-field effects [47j here and in the following). For 
a dipole at a distance z and oriented parallel to the 
mirror, we find from Eqs. © and HID 


Tde(z,w) = 7se,hM i 1 - 2 


sin(a) cos(a) sin(a) 


a 


a° 


(14a) 
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in terms of a = 2nuiz/c and the homogeneous-medium 
spontaneous-emission rate 7 S6i h = ii 2 noj 3 /( 37 rfe 0 c 3 ), i.e., 
n times the spontaneous-emission rate in free space. 
For a dipole emitter oriented normal to the mirror, the 
position-dependent spontaneous emission rate becomes 

| . (14b) 

In the limit z —> 0, the rate 7 ,^( 2 , u>) vanishes, while 
7 ^(z,w) tends to 27 sej h- In the limit 2 —> oo, both 
7 J e ( 2 ,u;) and 7 ,^( 2 , w) tend to the homogeneous-medium 
rate 7 se ,h(w). These two expressions are well known (hflj . 
Here we see how these exact results also follow from our 
multiple-scattering approach; we will compare their spa¬ 
tial dependence with that of the analytically determined 
energy transfer rates obtained by the same approach. 


liObw) = 7se,h(w) 0-3 


cos a 


sina 


III. ENERGY TRANSFER NEAR A MIRROR: 

PHENOMENOLOGY 

Figure [5] shows the total energy transfer rate between 
a donor and an acceptor as a function of distance 2 to the 
mirror. In this figure and all others below we use index 
n=l. The panels show results for several donor-acceptor 
spacings rd a = A/100, A/50, A/25. In all cases, the to¬ 
tal energy transfer reveals a considerable dependence at 
short range. In the limit of vanishing dipole-mirror dis¬ 
tance (z —> 0 ), dipoles perpendicular to the mirror have a 
four-fold enhanced transfer rate compared to free space. 
The factor four can be understood from the well-known 
method of image charges in electrodynamics: at a van¬ 
ishing distance, each image dipole enhance the field two¬ 
fold, and since energy transfer invokes two dipoles, the 
total result is a four-fold enhancement. 

With increasing dipole-mirror distance, the rate shows 
a minimum at a characteristic distance that is remarkably 
close to the donor-acceptor spacing 2 : ~ rda- At larger 
distances 2 > rd a , the transfer rate converges to the rate 
in the homogeneous medium. In Appendix [A] it is shown 
that this holds more generally: away from surfaces or 
other inhomogeneities, the FRET rate in an inhomoge¬ 
neous medium scales increasingly as oc l/(?r 4 r® a ), with n 
the refractive index surrounding the donor-acceptor pair. 
Figure [2] also shows that in the limit of vanishing dipole- 
mirror distance (2 —> 0 ), dipoles parallel to the mirror 
have an inhibited transfer rate. This result can also be 
understood from the method of image charges, since each 
image dipole provides complete destructive interference 
in the limit of zero distance to the mirror. 

With increasing dipole-mirror distance, the rate in¬ 
creases monotonously, and reaches half the free-space 
rate at a characteristic distance that is also remark¬ 
ably close to the donor-acceptor spacing 2 ~ rda- At 
larger distances 2 > ?'d a , the transfer rate tends to the 
homogeneous-medium rate. It is remarkable that even in 
a simple system studied here, a dramatic modification of 
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Figure 2: (Color online) Total energy transfer rate between a 
donor and an acceptor dipole, scaled to the free-space trans¬ 
fer rate, versus distance to the mirror, for the parallel and 
perpendicular configurations. The lower abscissa gives the 
distance in scaled units, and the upper abscissa absolute dis¬ 
tances at a wavelength A = (27r • 100)nm = 628nm. From 
top to bottom the three panels correspond to donor-acceptor 
spacings r<j a = A/100, A/50, A/25, where dipole-mirror dis¬ 
tances equal to ra a are marked by vertical dotted lines (off 
scale in the lowest panel). 


the energy transfer rate is feasible. In other words, we 
can already conclude that the energy transfer rate be¬ 
tween donor and acceptor is controlled by the distance 
to the mirror. The open question is whether this control 
is mediated by the LDOS. 

Figure |3] shows typical distances that characterize the 
distance dependent energy transfer rates in Figure [5] ver¬ 
sus donor-acceptor distance. For the perpendicular con¬ 
figuration we plot the distance where the transfer rate 
has a minimum, and for the parallel case we plot the dis¬ 
tance where the transfer rate equals 1 /2 of the free space 
rate. Both characteristic distances increase linearly with 
the donor-acceptor distance with near-unity slope. This 
behavior confirms that the distance dependence of the 
energy transfer rates in Figure [2] occurs on length scales 
comparable to the donor-acceptor distance. 
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Figure 3: (Color online) Characteristic mirror separations z/\ 
at which the total energy-transfer rate converges to the one in 
free-space. For dipoles perpendicular to the mirror, the char¬ 
acteristic distance is shown at which the total energy transfer 
rate 'yda/"/da,o has a minimum (see Fig. [2J. For dipoles par¬ 
allel to the mirror, the characteristic distance is shown where 
7da/7da,o equals 1/2 (see Fig. [2]). The lines are linear fits 
through the origin with slopes 3.75 and 3.45, respectively. 


Figure 0] shows the distance-dependence of the energy 
transfer rate in comparison to the spontaneous-emission 
rate. The latter varies with distance to the mirror on 
length scales comparable to the wavelength of light, as 
first discovered by Drexhage [U|. In contrast, the en¬ 
ergy transfer rates vary on dramatically shorter length 
scales, about one-and-a-lialf (parallel configuration) to 
two (perpendicular configuration) orders of magnitude 
smaller than the wavelength scale. 

To graphically investigate a possible relation between 
energy transfer rate and LDOS, Figure [5] shows a para¬ 
metric plot of the energy transfer rate as a function of 
(donor-only) spontaneous-emission rate, where each data 
point pertains to a certain distance z. The top abscissa is 
the relative LDOS at the donor emission frequency that 
equals the relative emission rate. The results at lower 
emission rate correspond mostly to the parallel dipole 
configurations in Figs. [2] and [4] whereas the results at 
higher emission rate correspond to mostly to the perpen¬ 
dicular configurations in these figures. For three donor- 
acceptor distances (ru a = A/100, A/50, A/25) Figure [5] 
shows that the energy transfer rate is independent of the 
emission rate and the LDOS over nearly the whole range, 
in agreement with conclusions of Refs. [13, HfJ 0, Ea]. 
The energy transfer decreases fast near the low emission 
rate edge and increases fast near the high emission rate 
edge, both of which correspond to distances very close 
to the mirror (cf. Fig. [2]). From Figure [5] it is readily 
apparent that the energy transfer rate does not increase 
linearly with the LDOS, leave alone quadratically, as pro¬ 
posed previously. The absence of a correlation between 
energy transfer rate and LDOS that is phenomenologi¬ 
cally shown here is one of our main results, and will be 


Figure 4: (Color online) Comparison of donor-acceptor en¬ 
ergy transfer rates 7d a and donor-only spontaneous emission 
rates 7 se , as a function of the distance 2 to the mirror. The 
lower abscissa is the scaled distance, the top abscissa is the 
absolute distance for A = 628 nm, both on a log scale. The 
energy transfer is scaled by the free-space energy transfer rate 
7da,o, the spontaneous emission by the free-space rate 7 se ,o- 
Data are shown both for the parallel and for the perpendicular 
configurations. For vanishing distance, 7d a /7da,o is inhibited 
to 0 for the parallel and enhanced to 4 for the perpendicular 
configuration. 


theoretically discussed in the remainder of this paper. 


IV. ENERGY TRANSFER VERSUS LDOS 


As is well-known, not all energy transfer is Forster en¬ 
ergy transfer. To find what fraction of the energy transfer 
corresponds to Forster energy transfer, we express the 
Green function in terms of the complete set of optical 
eigenmodes f A that satisfy the wave equation 

- V x V x f A (r) + E(r)(w\/ c) 2 f A (r) = 0, (15) 


with positive eigenfrequencies w A > 0. The Green func¬ 
tion, being the solution of Eq. ©, can be expanded in 
terms of these mode functions f A . An important prop- 
erty of this expansion can now be obtained from Ref. 4l| 
(in particular by combining Eqs. (21) and (22) of [4Tj). 
namely that the Green function can be written as the 
following sum of three terms: 


G(r,r» = c 2 ^ 


fA(r)f A *(rQ 
(w + ir/) 2 - uj 2 


(16) 
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Figure 5: (Color online) Parametric plots of scaled energy 
transfer rate versus scaled spontaneous-emission rate, for 
three donor-acceptor distances r da = A/100, A/50, A/25 from 
top to bottom, respectively. Data are obtained from Fig. [I] 
Black horizontal lines are constant at 7 da /7da,o = 1. The 
green dashed lines are linear relations 7 da = 7se. 


Since the Green function controls the energy transfer rate 
(see Eq. (0), it is relevant to discern energy transfer 
processes corresponding to these terms. The first term 
in Eq. (fl 6 l) denoted Gr corresponds to resonant dipole- 
dipole interaction (RDDI), the radiative process by which 
the donor at position r emits a field that is then received 
by the acceptor at position r'. The name ‘resonant’ de¬ 
scribes that photon energies close to the donor and ac¬ 
ceptor resonance energy are the most probable energy 
transporters, in line with the denominator (w + i ??) 2 — w 2 
of this first term. The second term in (fTTH) called Gg cor¬ 
responds to the static dipole-dipole interaction (SDDI) 
that also causes energy transfer from donor to acceptor, 
yet by virtual intermediate processes (see also Sec. V.B 
of Ref. fill) . As explained below, it is this SDDI that 
gives rise to the FRET rate that characteristically scales 
as r7 6 in homogeneous media and dominates the total en¬ 
ergy transfer for strongly subwavelength donor-acceptor 
separations. The third term in Eq. m is proportional to 


the Dirac delta function <5(r — r'). Since r ^ r' in case of 
energy transfer, this contribution vanishes. 

The fact that the Green function can be written as 
the sum of three terms as in Eq. m is important, and 
implies that for arbitrary environments the static part of 
the Green function can be obtained from the total Green 
function by the following limiting procedure (for r ^ r') 

Gs(r,r',w) = \ lim w 2 G(r, r', w), (17) 

to w-m 

which provides a justification of our use of the term 
‘static’. As an important test, selecting in this way 
the static part of the Green function of a homogeneous 
medium 0 indeed gives that only 

G hlS (ri,r 2 ,u;) = — 5^3 0 - 3fr), (18) 

with r = r d — r 2 contributes to Forster energy trans¬ 
fer, and not the terms of Gh that vary as 1 /r and 1 /r 2 . 
Incidentally and by contrast, for general inhomogeneous 
media the static Green function does not only depend on 
the distance between donor and emitter, but rather on 
the absolute positions of both donor and acceptor in the 
medium. 

Having thus defined Forster energy transfer as that 
part of the total energy transfer that is mediated by the 
static dipole-dipole interaction, we can now also define 
the square of the Forster transfer amplitude, in analogy 
to Eq. 0, by 


w F ( r a ,r d ,u;) = ^ ' G s (r a ,r d ,w) • fj, d \ 2 . 

(19) 

This appears similar to Eq. m, yet with the total Green 
function G replaced by its static part Gg, as defined 
in Eq. (fTbl) and computed in Eq. (fI71) . The FRET 
rate 7f is then obtained by substituting WF(r a , r d , w) for 
w(r a ,r d ,w) into Eq. 0, giving: 

/*00 

7 F (r a ,r d )= / dw<T a (w)w F (i , a,rd,w)<T d (w). (20) 


Here we arrive at an important simplification in the de¬ 
scription of Forster transfer in inhomogeneous media, 
by noting that from Eqs. o and (HU), the quantity 
WF(ra, I'd, w) is independent of frequency ui. The FRET 
rate 7 f is then given by the simple relation 

/ OO 

dwCTaM CT d (w). (21) 

-OO 


While this expression looks similar to the approximate 
expression for the total energy transfer rate (Eq. []}, we 
emphasize as a first point that Eq. m is an exact expres¬ 
sion for the FRET rate, even for broad donor and accep¬ 
tor spectra. A second crucial point is that the spectral 
overlap integral in Eq. m is the same for any nanopho¬ 
tonic environment [60}. All effects of the nondispersive 
















inhomogeneous environment are therefore contained in 
the frequency-independent prefactor wp(r a , r d ). In other 
words, while there is an effect of the nanophotonic envi¬ 
ronment on the FRET rate (see Fig. n, this effect de¬ 
pends only on the donor and acceptor positions but does 
not depend on the resonance frequencies of the donor and 
acceptor (for constant medium-independent overlap inte¬ 
gral in Eq. ((21)1 ). If the FRET rate does not depend on 
the donor and acceptor frequencies, then the FRET rate 
can not be a function of the LDOS at these particular 
frequencies. A third crucial point is that this conclusion 
is valid for any photonic environment that is lossless and 
weakly dispersive in the frequency range where the donor 
and acceptor spectra overlap, hence this conclusion is not 
limited to an ideal mirror. 

For homogeneous media it is well known that Forster 
energy transfer dominates the total energy transfer at 
strongly sub-wavelength distances, and we will now see 
that this is also the case in inhomogeneous media, again 
taking the ideal mirror as an example. The total energy 
transfer near an ideal mirror depends on the total Green 
function as given in Eqs. m and (fl3l) for the two dipole 
configurations (cf. Fig. [T]) . For the donor and acceptor 
near the mirror in the parallel configuration, we use the 
procedure of Eq. G3 and obtain for the static parts 


/x ll -G s (r a ,r d ,w)-/x li = 


H 2 c 2 


4 \ r(j a (yfi^+Zz 3 ) 3 j ' 

( 22 ) 

while for the perpendicular configuration we find 


^ ■ G s (r a ,r d ,w) • = 


M 2 c 2 

4-7m 2 u; 2 1 r da 


(VrL + 4* 2 ) 3 


1-3- 


4.z 2 


da 


(23) 


We note that in both cases the static interaction in a ho¬ 
mogeneous medium is recovered for FRET pairs at dis¬ 
tances to the mirror much larger than the donor-acceptor 
distance z r da . This agrees with the large-distance 
limits shown in Figs. [2] and |U The spatial dependence 
of the Forster transfer amplitude of Eq. m and of the 
FRET rate in Eq. (1211) is hereby determined for both 
conhgur ations. 

In Figure (HI we display the ratio of the FRET rate 
and the total energy transfer rate as a function of donor- 
acceptor distance, for three distances z of the FRET pair 
to the mirror, and for both dipole configurations, taking 
n=l. For the total rate we use the narrow bandwidth 
assumption of Eq. a. Irrespective of the distance to the 
mirror and of dipole configuration, the total energy trans¬ 
fer rate equals the FRET rate for all practical purposes, 
as long as (r da w/c -C 1), in other words for strongly sub¬ 
wavelength donor-acceptor distances. Intriguingly, with 
increasing donor-acceptor distance beyond typical FRET 
distances, the ratio of the two rates exceeds unity. Since 
the total rate 7 da equals the absolute square of the sum of 
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Figure 6: (Color online) Forster resonance energy transfer 
rate scaled to the total energy transfer rate (7F/7 da ) versus 
donor-acceptor distance r da for three distances 2 of donor and 
acceptor to the mirror. The upper panel is for dipoles par¬ 
allel to the mirror, the lower panel for dipoles perpendicular 
to the mirror. Note the logarithmic r da , with dimensionless 
scaled values on the lower abscissa and absolute distance in 
nanometers on the upper abscissa for A = 628nm. 


static and resonant transfer amplitudes, and the two am¬ 
plitudes interfere destructively in this intermediate dis¬ 
tance range, the FRET (static) rate can indeed dominate 
the total energy transfer rate. At large donor-acceptor 
distances (r da w/c 1), the FRET rate decreases much 
faster with distance than the total transfer rate, similar 
as in homogeneous media. In this large-distance range, 
the energy transfer is radiative: the donor emits a photon 
that is absorbed by the acceptor. 

Figure [7] is complementary to the previous one in the 
sense that here the FRET rate is plotted versus dis¬ 
tance to the mirror 2 for several donor-acceptor distances 
r da , and for both dipole configurations. We again show 
the ratio of the FRET rate and the total transfer rate 
(using Eq. (HJ) for 7 da ). At donor-acceptor distances 
r da = A/100 and r da = A/50, typical for experimental 
situations, we clearly see that FRET dominates the total 
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Figure 7: (Color online) FRET rate 7 f divided by the total 
energy transfer rate 7 da , versus distance to the mirror, for 
three values of the donor-acceptor distance r da . The lower 
abscissa is the dimensionless reduced distance, the upper ab¬ 
scissa is the absolute distance in nanometer for A = 628nm. 
Upper panel: parallel dipole configuration; lower panel: per¬ 
pendicular dipole configuration. 


energy transfer rate, independent of the distance to the 
mirror. At least 98% of the total energy transfer rate con¬ 
sists of the FRET rate. Even for a large donor-acceptor 
distance r da = A/20 that is much larger than in most 
experimental cases (that is, r da = 31nm at A = 628 nm), 
the FRET rate and the total rate differ by only some 
ten percent. Thus, Figures [6] and [7] illustrate that in 
the nanophotonic case near an ideal mirror, the FRET 
greatly dominates the total energy transfer at strongly 
sub-wavelength donor-acceptor distances, similar as in 
the well-known case of homogeneous media. Therefore, 
we conclude that not only is there no correlation between 
the LDOS and the total transfer rate, there is also no cor¬ 
relation between the LDOS and the FRET rate either. 


V. ENERGY TRANSFER IN TERMS OF A 
FREQUENCY-INTEGRATED LDOS 

Our general derivation of the FRET rate in a weakly 
dispersive nanophotonic medium (Eq. 1211) has convinc¬ 
ingly shown that the FRET rate has no dependence on 
the local density of optical states evaluated at the donor’s 
resonance frequency. In this section we will insist on es¬ 
tablishing a link between the FRET rate and the LDOS, 
if only to counter the argument that we are from the 
outset biased against such a relation. Interestingly, the 
relation that we derive will also shed a new light on efforts 
to control the FRET rate by LDOS engineering. 

We start with the mode expansion of the Green func¬ 
tion in Eq. ITCl) to derive a useful new expression, relat¬ 
ing the Forster transfer rate to a frequency-integral over 
Im[G], We use the fact that Gs(r, r',w) is real-valued, 
as is proven in Ref. [49] . Thus the imaginary part of the 
Green function is equal to Im[Gu] and the mode expan¬ 
sion of Im[G] becomes 

Im[G(r,r',w)] = ^ f A (r)f^(r')<5(a; - w A ), (24) 

A 

with to > 0. We note that only modes with frequencies 
uj\ = lu show up in this mode expansion of Im[G], This 
can also be seen in another way: the defining equation 
for the Green function Eq. © implies that the imaginary 
part of the Green function satisfies the same source-free 
equation m as the subset of modes f A (r) for which the 
eigenfrequency oj\ equals oj. Therefore, Im[G(r, r',w)] 
can be completely expanded in terms of only those de¬ 
generate eigenmodes. The mode expansion (l24ll is indeed 
a solution of Eq. (flbl) . If we multiply the right-hand side 
of Eq. (1241) with u> and then integrate over w, we obtain 
as one of our major results an exact identity for the static 
Green function Gs 

2 

G s (r a ,r d ,w) = —y/ dwiWiIm[G(r a ,r d ,wi)]. (25) 

^ Jo 

This identity is valid for a general nanophotonic medium 
in which material dispersion can be neglected. Eq. © 
was derived using a complete set of modes, yet does not 
depend on the specific set of modes used. When inserting 
this identity into Eq. (HI), we obtain the Forster transfer 
rate u>s(oj) and hence the transfer rate 7 f of Eq. m 
in terms of the imaginary part of the Green function. 
While this is somewhat analogous to the well-known ex¬ 
pression for the spontaneous-emission rate (Eq. [5|), there 
are two important differences: The first difference be¬ 
tween Eq. (1251) for Forster energy transfer and Eq. © 
for spontaneous emission in terms of Im[G] is of course 
that Eq. (1251) is an integral over all positive frequencies. 
The second main difference is that in Eq. C5l) the Green 
function Im[G(r a , r d , wi)] appears with two position ar¬ 
guments - one for the donor and one for the acceptor - 
instead of only one position as in the spontaneous emis¬ 
sion rate. A major advantage of an expression in terms 
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of Im[G] is that Im[G] does not diverge for r a —> r d , in 
contrast to Rc[G], 

In Appendix |B] we verify and show explicitly that the 
identity in Eq. (1251) holds both in homogeneous media 
as well as for the nanophotonic case of arbitrary po¬ 
sitions near an ideal mirror. The identity in Eq. (1251) 
is important since it is our goal in this section to ex¬ 
plore whether FRET rates are functionally related to the 
LDOS in any way. With Eq. (ESI), both quantities can 
now be expressed in terms of the imaginary part of the 
Green function, which brings us considerably closer to 
the goal. 

We now use Eq. (l25l) to derive an approximate expres¬ 
sion Gg L * for the static Green function Gg that allows 
us to relate the Forster transfer rate to the frequency- 
integrated LDOS. Our approximation is motivated by 
the fact that Im[G(rd — r a ,w)] for homogeneous media 
(based on Eq. (flrjl) varies appreciably only for varia¬ 
tions in the donor-acceptor distance r da on the scale of 
the wavelength of light, typically rd a ~ Ao = 500 nrn 
(with Ao = 2ttc/loo). From Eq. (fill) it follows that the 
same holds true for Im[G(r d , r a , w)] for the ideal mir¬ 
ror. In contrast, Forster energy transfer occurs on a 
length scale of rd a — 5 nm, typically a hundred times 
smaller. Motivated by these considerations, we approx¬ 
imate Im[G(r a ,r d ,uq)] in the integrand of Eq. (E5l) by 
the zeroth-order Taylor approximation Im[G(rd, rd, uq)]. 
The accuracy of this approximation depends on the opti¬ 
cal frequency ui. The approximation will not hold for all 
frequencies that are integrated over, and becomes worse 
for higher frequencies. But it appears that we can make 
an accurate approximation throughout a huge optical 
bandwidth 0 < uq < O. If we choose O = 10wo, be., 
a frequency bandwidth all the way up to the vacuum 
ultraviolet (VUV), then Im[G(r a , r d , oq)] will only devi¬ 
ate appreciably from Im[G(r d , r d , uq)] for donor-acceptor 
distances rd a > Ao/10, which is in practice of the order of 
50 nm, much larger than typical donor-acceptor distances 
in Forster transfer experiments. We obtain the expres¬ 
sion for the approximate static Green function Gg L ^ as 

Gg L) (r a ,r d ,w) = — - / duq uq Im[G(r d , r d , uq)] 

-xu 1 J 0 

2 

+ - 2 / duq uq Im[G(r a , r d , uq)(]26) 

7ruJ J n 

Here the first term is recognized to be an integral of 
the LDOS over a large frequency bandwidth, ranging 
from zero frequency (or ‘DC’) to high frequencies in 
the VUV range. While the specific value of D does 
not matter much, it is important that D can be chosen 
much greater than optical frequencies, while the inequal¬ 
ity n(r d )r2r da /c <C 1 still holds. Within this approxi¬ 
mation, the Forster transfer rate can be related to the 
frequency-integrated LDOS: if we replace Gg by Gg L ^ in 
Eq. (fl9l) for Wp , thereby obtaining the LDOS approxi¬ 
mation Wp L * (r a , r d , u) = |/i* ■ Gg L (r a ,r d ,u) • /i d | 2 , and 
approximate wf in Eq. (1201) by this , we obtain the 
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Figure 8: (Color online) LDOS-approximated FRET rate 7 p L) 
(Eq. (1271) ) normalized to the exact FRET rate 7 f (Eq. (1201) ) 
versus the bandwidth D of the LDOS-frequency integral. 
Lower abscissa: 0 scaled by the donor frequency u d = 2nc/X. 
Upper abscissa: minimum wavelength A m i n = 2-kc/^I for 
A = 628nm. Black full curves are for dipole-to-mirror distance 
z = A/100, red dashed curves for z = A/40, and blue dashed- 
dotted curves for a = A/2, all curves are for a donor-acceptor 
distance r da = A/100, (a) Parallel dipole configuration; (b) 
perpendicular dipole configuration, (c) Comparison of the 
LDOS-approximation (Eq. (1271) ) and the high-frequency ap¬ 
proximation (Eq. (1281) ) of the FRET rate as a function of 
LDOS bandwidth Q. Rates are scaled to the exact FRET 
rate, and the distance to the mirror and the donor-acceptor 
distance are fixed. 


relation between the approximate FRET rate 7 ^ and 
the LDOS to be [U] 

/ OO 

dwCT a (w) w^ L) (w)a d (w). (27) 

-OO 

In Figure [5] we verify the accuracy of the LDOS- 
approximated FRET rate 7 ^ near the ideal mirror, by 
varying the frequency bandwidth Q over which we make 
the approximation. The required frequency integrals of 
Eq. (E^l) are calculated analytically in Appendix 1C 21 In 
Fig. [5] we see that for both dipole configurations, the ap- 
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Figure 9: (Color online) LDOS-approximated FRET rate 
(Eq. (1271) 1 normalized to the exact FRET rate 7 f (Eq. (1201) 1 
versus (scaled) distance to the mirror, for an LDOS frequency 
bandwidth up to ft = lOwo (red dashed curve), and up to 
fl = 2cdo (blue dashed-dotted curve). The donor-acceptor 
distance is r da = A/100. 


proximate FRET rate indeed tends to the exact rate for 
vanishing fl. For fl up to 10w d , the approximate rate is 
very close to the exact one, to within 5%. Even at higher 
frequencies, up to fl = 20a7 d , the approximate FRET 
rate is within 10% of the exact rate, as anticipated on 
the basis of our general considerations above. 

The validity of the approximate FRET rate 7^ im¬ 
proves when the donor-acceptor distance r da is reduced, 
since the spatial zero-order Taylor expansion of Im[G] is 
then a better approximation. We can also improve the 
approximation by reducing the frequency bandwidth fl in 
which we make the Taylor approximation. Both trends 
are indeed found in Anoendix IC ll where 7^ is calculated 
for the homogeneous medium. In the limit of a vanishing 
frequency bandwidth (fl —>- 0), the approximate Forster 
transfer rate jp reduces to the exact Forster transfer 
rate 7 f of Eq. (l20l) . 

To verify that the approximate FRET rates shown in 
Fig.[S]were not ‘lucky shots’ for the chosen fixed distances 
to the mirror, we study in the complementary Figure[5|a) 
the accuracy of 7^ as a function of distance to the mir¬ 
ror z , for a constant LDOS bandwidth fl = lOwd- The 
figure clearly shows the great accuracy of the LDOS ap¬ 
proximation, irrespective of the distance z of the FRET 
pair to the ideal mirror. For a narrower bandwidth of 
fl = 2w d , the accuracy is even better, as expected. 

At this point, one may naively conclude from Figures[8] 
and [9] that the FRET rate is intimately related to an in¬ 
tegral over the LDOS. This is too rash, however, because 
the corresponding approximate relation Eq. (ED) consists 
of two integrals, where only one of them is an integral 
over the LDOS, while the other is a high-frequency in¬ 
tegral of the complex part of the Green function. Thus 
the relevant question becomes: what happens if we make 


a cruder approximation to the FRET rate by simply re¬ 
moving the LDOS integral? Instead of Eq. (l26|) we then 
use the high-frequency approximation (HF) to the static 
Green function Gg HF *: 

Gf F) (r a ,rd,w) = ^7 / dwi uq Im[G(r a , r d ,aq)]. 
nuj ~ J n 

(28) 

This leads to a high-frequency approximation 

/TT-p\ 

for the squared Forster amplitude Wp (lo) = 
{2Tr/h 2 ){u]/(E 0 c 2 )) 2 \nl ■ G s HF) ( r a: r d,w) • /x d | 2 , and 
a high-frequency approximation the FRET rate 7p ; : 

/ OO 

dwcr a (w)«;|, HF) (r a ,r d ,a;)cr d (a;). (29) 

-CO 

In Figure (Hfc) the two approximated FRET rates 7^ 

and 7p ; are compared for the ideal mirror, both scaled 
by the exact FRET rate yp, as a function of the band¬ 
width fl. The donor-acceptor distance and the distance 
to the mirror are fixed. Indeed 7^ is the more accurate 
approximation of the two, yet 7p HF ^ is not a bad ap¬ 
proximation at all: by only integrating in Eq. (1281) over 
high frequencies or > fl = 10w d , 7p ; is accurate to 
within about 7%. If we take a narrower - yet still broad 
- frequency bandwidth, for example up to fl = 2w d (in 
the UV), we still neglect the LDOS in the whole visible 
range. Nevertheless Figure [51c) shows that for fl = 2w d 
the two approximations 7^ and 7p HF ' 1 agree to a high ac¬ 
curacy with the exact rate yp. Therefore, Figures [8] and 
[9] show that for the ideal mirror there is essentially no de¬ 
pendence of the FRET rate on the frequency-integrated 
LDOS at visible frequencies, and only a weak dependence 
on the frequency-integrated LDOS at UV frequencies and 
beyond. We note that this conclusion is different from 
and complementary to the one in Sec. El where the 
FRET rate was found to not depend on the LDOS at 
one frequency, namely at the transferred energy fiw d . 


VI. DISCUSSION 

In this section, we discuss consequences of our theoreti¬ 
cal results to experiments, first regarding relevant length 
scales. We have performed analytical calculations and 
plotted rates versus reduced lengths, namely the reduced 
distance to the mirror zuj/c = 2nz/\ , and the reduced 
donor-acceptor distance r da /A. To increase the relevance 
of our results to experiments and applications, we have 
plotted in several figures additional abscissae for abso¬ 
lute length scales that pertain to a particular choice of 
the donor emission wavelength A d - Here we have chosen 
A d = 2-7T /u)d = (27r • 100)nm ~ 628 nm, a figure that we 
refer to as a ’’Mermin-wavelength” [53j, as it simplifies 
the conversion between reduced units and real units to a 
mere multiplication by 100 x. 
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Figures [2] and [3] characterize the distance dependence 
to the mirror. The range where both the total and the 
FRET rate are controlled by the distance to the mirror is 
in the range z < 20 nm. This range is set by the donor- 
acceptor distance that is for most typical FRET pairs in 
the order of r& a = lOnm, in view of typical Forster dis¬ 
tances of the same magnitude [1]. Interestingly, while the 
energy transfer is in this range (z < 7'da) not controlled 
by the LDOS, the transfer rate itself is nevertheless con¬ 
trolled by precise positioning near a mirror. An example 
of a method that could be used to achieve such control 
at optical wavelengths is by attaching emitters - such as 
molecules or quantum dots - to the ends of brush poly¬ 
mers with sub-10 nm lengths [HI]. With Rydberg atoms, 
it appears to be feasible to realize the situation z < rda, 
albeit in the GHz frequency range [55| . 

Figure [6] characterizes the donor-acceptor distance de¬ 
pendence of the transfer rate. It is apparent that Forster 
transfer dominates in the range rd a < 20nm, a length 
scale much smaller than the wavelength of light. In the 
range rda > 100 nm, energy transfer is dominated by 
radiative transfer, which is reasonable as this distance 
range becomes of the order of the wavelength. 

Let us now briefly discuss broadband LDOS control. If 
one insists on invoking the LDOS to control the FRET 
rate, Figure [5] shows that one must control the LDOS 
over a huge bandwidth that ranges all the way from zero 
frequency (’DC’) to a frequency D that is on the order of 
10 times the donor emission frequency uid- This agrees 
with the qualitative statements in Ref. [2?J. If we con¬ 
sider the Mermin-wavelength 628 nm, the upper bound 
on the LDOS bandwidth corresponds to a wavelength 
of 62.8nm nm, which is deep in the vacuum ultraviolet 
(VUV) range. At these very short wavelengths, all mate¬ 
rials that are commonly used in nanophotonic control - be 
it dielectrics such as silica, semiconductors such as silica, 
or metals such as silver or gold - are strongly absorbing. 
In practice, the optical properties of common nanopho¬ 
tonic materials deviate from their commonly used prop¬ 
erties at wavelengths below 200 to 250 nm, which cor¬ 
responds to O < 3 uid- Yet, even if one were able to 
control the LDOS over a phenomenally broad bandwidth 
0 < D < 3uJd, Figure [9] shows that the broadband LDOS- 
integral contributes negligibly - much less than 10“ 3 - to 
the Forster transfer rate. Thus, with the current state of 
the art in nanofabrication, true LDOS-control of Forster 
energy transfer seems to be extremely challenging. 

The importance of distinguishing FRET from other 
energy-transfer mechanisms has also been emphasized 
by Govorov and co-workers [51 j , who studied plasmon- 
enhanced Forster transfer near conducting surfaces. 
They predict that near metal surfaces the FRET rate 
can be enhanced by much more than the factor 4 reported 
here for an ideal mirror, but that a strong enhancement 
occurs only near the plasmon peak, i.e., in a highly dis¬ 
persive region, while on the other hand not too close to 
the resonance since otherwise loss makes Forster trans¬ 
fer invisible. Analogous FRET enhancements only in the 


highly dispersive region near a resonance were found in 
Ref. [25j ] for single-resonance Drude-Lorentz type media. 
In Ref. [52j, FRET near graphene is also clearly distin¬ 
guished from long-range plasmon-assisted energy trans¬ 
fer. As a future extension of our work, it will be inter¬ 
esting to study plasmon-enhanced FRET and long-range 
energy transfer rates by considering a mirror with a res¬ 
onance. 

How do our theoretical results compare to experi¬ 
ments? Our theoretical findings support the FRET-rate 
and spontaneous-emission rate measurements by Blum et 
al. [30], where it was found that Forster transfer rates are 
unaffected by the LDOS. Our findings also agree with the 
results of Refs. [13, HI, [H, Hsl • What about other exper¬ 
imental studies that do report a relation between FRET 
rate and LDOS? Assuming our theory to be correct, this 
discrepancy can mean three things. First , it could be 
that in those experiments the energy transfer between 
donor and acceptor separated by a few nanometers was 
not dominated by Forster transfer. This might occur 
for energy transfer within a high-Q cavity, but otherwise 
does not seem to be a probable explanation. Second , our 
theory leaves open the possibility that there is a corre¬ 
lation between FRET rates and LDOS in case of strong 
dispersion and/or loss in the frequency overlap range of 
donor and acceptor spectra, in case of plasmon-mediated 
FRET for example, because in that case our theory does 
not apply. However, if the correlation between LDOS 
and FRET rates indeed relies on strong dispersion, then 
this correlation would be a particular relation rather than 
the sought general relation. Moreover, Ref. HH theoreti¬ 
cally studies spontaneous-emission and FRET rates near 
metal surfaces but does not report a general linear rela¬ 
tion between them. Technically, in case of non-negligible 
absorption and concomitant complex dielectric function, 
the concept of the LDOS breaks down, but FRET rates 
can still be compared to the imaginary part of the Green 
function. As a third possible reason why our theory does 
not predict a relation between FRET rates and LDOS 
while some experiments do, we should mention that our 
theory does not include typical aspects of experiments, 
such as incompletely paired donors, cross-talk between 
dense donor-acceptor pairs, inhomogeneously distributed 
donor-acceptor distances, or transfer rates influenced by 
strongly inhomogeneous field distributions that may oc¬ 
cur near nanoparticles or slits and indentations. While 
describing these effects adds substantial complexity to 
our simple model, they may be taken into account by 
judicious choices of simple scatterers such as spheres or 
dipoles. Nevertheless, it does not seem likely that such 
particular additional effects will induce a general depen¬ 
dence of the FRET rate on the LDOS. 

Regarding the subject of quantum information process¬ 
ing, FRET is a mechanism by which nearby (< 10 nm) 
qubits may interact p~5l-f20| . intended or not. Lovett et 
al. E3 considered the implications of Forster resonance 
energy transfer between two quantum dots. In one im¬ 
plementation, it was found that it is desirable to sup- 
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press the Forster interaction in order to create entangle¬ 
ment using biexcitons. In another implementation, it was 
found that Forster resonance energy transfer should not 
be suppressed, but rather switched. There is a growing 
interest in manipulating the LDOS, either suppressing it 
by means of a complete 3D photonic band gap 1331 . or 
by ultrafast switching in the time-domain [561 ]. It fol¬ 
lows from our present results that these tools cannot be 
used to also switch or suppress Forster resonance energy 
transfer between quantum bits. On the other hand, our 
results do indicate that FRET-related quantum informa¬ 
tion processing may be controlled by carefully positioning 
the interacting quantum systems (*.e., the quantum dots) 
in engineered inhomogeneous dielectric environments. 


VII. CONCLUSIONS 

Using an exactly solvable analytical model, we have 
seen that Forster resonance energy transfer rate from a 
donor to an acceptor differs from the one in a homoge¬ 
neous medium in close vicinity of an ideal mirror. For 
two particular dipole configurations, we found that the 
FRET transfer rate is inhibited (to 0) or markedly en¬ 
hanced (4x). Thus, even this simple model system of¬ 
fers the opportunity to control energy transfer rates. It 
turns out that differences in FRET rates as compared 
to those in a homogeneous medium are only noticeable 
at distances to the mirror on the order of the donor- 
acceptor distance rd a or smaller. This distance together 
with the wavelength are the only natural length scales 
in this simple problem. Since r da is typically less than 
10 nm, that is, orders of magnitude smaller than an op¬ 
tical wavelength, any substantial variations in the FRET 
rates due to the nanophotonic environment occur on a 
distance scale on which the LDOS does not vary appre¬ 
ciably in the dielectric medium near the ideal mirrror. 

On the larger distance scale of an optical wavelength 
away from the mirror, there are well-known LDOS vari¬ 
ations. At these larger distances, the FRET rate is con¬ 
stant and the same as in a homogeneous medium. So the 
one quantity varies appreciably when the other does not, 
and vice versa. Therefore, we conclude that the FRET 
rate does not correlate with the partial or total LDOS. 
This particular example already suffices to conclude that 
in general , the FRET rate does not correlate - neither lin¬ 
early, nor quadratically, or otherwise - with the LDOS. 

How large is the class of environments for which FRET 
rate and LDOS are uncorrelated? We have derived as 
one of our main results the simple and exact expres¬ 
sion (Eq. l2ll) for the FRET rate in an inhomogeneous 
medium. As a consequence, it follows that Forster en¬ 
ergy transfer rates are independent of the LDOS at the 
transferred photon frequency in all nanophotonic media 
where material dispersion and loss can be neglected in the 
donor-acceptor frequency overlap interval. Other main 
results are the exact expression (Eq. [55]) of the static 
Green function in terms of a frequency integral over the 


imaginary part of the total Green function, and the corre¬ 
sponding approximate relation (Eq. [271) between FRET 
rate and the frequency-integrated LDOS. We used the 
latter relation to show that FRET rates near an ideal 
mirror are numerically independent of the LDOS, even 
when integrating the LDOS over all visible frequencies. 
We have also argued why the same will be true for other 
media with weak material dispersion as well. 

We have emphasized that not all energy transfer is 
Forster energy transfer, and that for typical Forster 
donor-acceptor distances below 10 nm, the energy trans¬ 
fer is dominated by Forster transfer, as the example of 
the ideal mirror has also shown. For arbitrary photonic 
environments we defined Forster transfer as being medi¬ 
ated by virtual photon exchange, the strength of which 
is determined by the static Green tensor, which in homo¬ 
geneous media gives rise to the characteristic 1 / (n 4 r® a )- 
dependence of the Forster resonance energy transfer rate. 
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Appendix A: Scaling with donor-acceptor distance 
of Forster transfer rate 

Here we show that the homogeneous-medium Forster 
transfer rate, scaling as oc l/(nj 1 1 7‘® a ), is an important 
limiting case also for inhomogeneous media. Let us as¬ 
sume that the donor and acceptor are separated by a 
few nanometers, experiencing the same dielectric mate¬ 
rial with a dielectric function £h> within an inhomoge¬ 
neous nanophotonic environment. In all of space, we de¬ 
fine the optical potential V(r, oj) = — [e(r) — £h](w/c) 2 l, 
so that the optical potential vanishes in the vicinity of 
the donor-acceptor pair. Then the Green function of the 
medium can be expressed in terms of the homogeneous- 
medium Green function and the optical potential as 

G(r a ,r d ,w) = G h (r a - r d ,w) (Al) 

+ /dn G h (r a - ri,w) • V(ri,w) • G(ri,r d ,w), 

which is the Dyson-Schwinger equation for the Green 
function that controls the energy transfer. The equa¬ 
tion can be formally solved in terms of the T-matrix of 
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the medium as 

G(r a , r d ) = G h (r a - r d ) (A2) 

+ j dridr 2 G h (r a - ri) • T(ri,r 2 ) • G h (r 2 - r D ), 

where the frequency dependence was dropped for read¬ 
ability. The important property of the T-matrix 
T(ri,r 2 ,w) is now that it is only non-vanishing where 
both V(ri) and V(r 2 ) are nonzero, so that it vanishes in 
the vicinity of the donor-acceptor pair. Thus the Green 
function that controls the energy transfer is given by the 
sum of a homogeneous-medium Green function and a 
scattering term. The former is a function of the distance 
between donor and acceptor, whereas the latter does not 
depend on the D-A distance, but rather on the distance 
of donor and acceptor to points in space where the optical 
potential is non-vanishing. 

As the donor-acceptor distance r da is decreased, the 
homogeneous-medium contribution in Eq. (IA2I) grows 
rapidly, essentially becoming equal to Gh,s(r a — r d ,w) 
of Eq. (USD, whereas the contribution of the scattering 
term does not change much. So in the limit of very small 
r da , or when making the distance to interfaces larger, the 
homogeneous-medium term always wins, and one would 
find the well-known Forster transfer rate of the infinite 
homogeneous medium oc l/(«h r da)- 

Appendix B: Tests of identity (1251) 

1. Test for a homogeneous medium 


J 0 °° dxx sin(fca;) = while the third term in 

the square brackets of Eq. (IB2I) does contribute since 
f 0 da:sin(a:)/a; = w/2. Thus the projection of the iden¬ 
tity mi indeed holds for spatially homogeneous media. 

2. Test for an ideal mirror 

For the ideal mirror we again only consider a projec¬ 
tion of the identity (l25l) . first projecting onto dipoles 
corresponding to the parallel configuration of Fig. [I] 
The Green function for the ideal mirror is given in 
Eq. (fill) , and its static part for the parallel configuration 
by Eq. (l22l) . Now for this parallel configuration, the pro¬ 
jected Green tensor consists of a homogeneous-medium 
and a reflected part, and so does the projected static 
Green function. In Sec. IB ll above we already showed that 
the sought identity indeed holds for homogeneous media. 
So the remaining task is to show that the identity C5l) 
holds separately for the reflected parts of the projected 
Green functions. This is not difficult since mathemati¬ 
cally the frequency integral that is to be performed is the 
same as for the homogeneous medium; only the distance 
parameter r da is to be replaced by \J r^ a + 4 z 2 . Thus the 
projection of the identity (E5l) onto the parallel dipole 
directions indeed holds. The qualitative novelty as com¬ 
pared to the homogeneous-medium case is that we thus 
show that the identity holds irrespective of the distance 
z of the FRET pair to the mirror. We also checked (not 
shown) that the identity C5l) holds for the projection 
onto perpendicular dipoles, i.e. as in the perpendicular 
configuration of Fig. [TJ 


The Green function Gh(r, ui) for homogeneous media 
is given in Eq. ©, and its static part by Eq. (fT51) . The 
identity (1251) that relates them can be shown to hold as 
a tensorial identity; here we derive the identity for its 
projection fi ■ Gh(r, ui) ■ fi, where we assume ji to be per¬ 
pendicular to r. (Physically, this corresponds to energy 
transfer between donor and acceptor with equal dipoles 
both pointing perpendicular to their position difference 
vector.) The projection of the identity (l?5l) that we are 
to derive has the form 


/r 2 c 2 


47m 2 w 2 r da 


27 T 2 n 2 L 0 2 r da 


dwiuqlm 


(Bl) 


e in “ rd “/ c P(mwr da /c) 


Now by integration variable transformation the right- 
hand side of this equation can be worked out to give 


Appendix C: Accuracy of the approximate 
expressions (1261) and (1291) 

To test the accuracy of the LDOS approximation 
Gg L ^(r a ,r d ,w) of the static Green function, it is conve¬ 
nient to use Eq. (E5l) to rewrite Eq. (E(jl) as 

G^ L) (r a ,r d ,w) = G s (r a ,r d ,w) (Cl) 

2 f n 

4-y / dwi wiIm[G(r d ,r d ,wi) - G(r a ,r d ,wi)], 

KUJ Jo 

In this form, the approximate static Green function is 
equal to the exact expression plus an integral over a fi¬ 
nite interval of a well-behaved integrand. Likewise, to 
test the accuracy of the high-frequency approximation 
Gg HF ^(r a ,r d ,w) defined in Eq. (l29l) of the static Green 
function, it is useful to rewrite it as 


H 2 c 2 


« w r da Jo 


,oo 

. . ,, , sm(x) 

cos(kx) + xsm(kx) - 

/ dx 

Jo 

X J 
( 


(B2) 


with dummy variable k equal to unity. Now the 
first two terms within the square brackets do not con¬ 
tribute to the integral since dxcos(kx) = nS(k ) and 


Gg HF) (r a ,r d ,w) = G s (r a ,r d ,w) (C2) 

2 f n 

- —5- / dwiWiIm[G(r a ,r d ,wi)], 
7TW JO 

Again the integrand is well-defined, i.e. non-diverging, 
over the entire finite integration interval. 
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1. Accuracy of LDOS approximation for 
homogeneous media 


2. Accuracy of LDOS approximation for the ideal 
mirror 


We estimate the accuracy of Eq. ei) for the Green 
function (JHl) of a homogeneous medium. By taking the 
projection onto dipole vectors both on the left and right, 
we find 


/*f2 

Im / dwiwi/x • [G h (r d ,r d ,wi) 
J o 


p 2 n 

47TC 



dujiujl 



Gh(r a ,r d ,wi)] • ft = 
(C3) 


where D = ntod/c and for convenience we defined the 
function h(x) = cos(x)/x 2 + sin(a:)(l — l/x 2 )/x. So here 
(and also for the mirror below) we must determine inte¬ 
grals of the type 


= / duiujlh(ujia) (C4) 

Jo 

= {Ft/A) 3 [sin(A) — Acos(A) — S't(A)]. 

where A = f la and Si [a:] = fj dtsin(t)/t is the sine inte¬ 
gral. For f la -C 1 we find the approximation 

"<“■») = ! fiS - < C5 > 

With this result, we find that the relative error of making 
the LDOS approximation G^g^a, r d , to) of Eq. (IC1I) for 
the Green function Gh i s(r a , r d , w) is 


M • ( G hi - G h,s) • ft 4 

- 7Tc - 7, - = • ( C6 ) 

ft • G h ,s • ft 757T 


This fifth-power dependence shows that for homogeneous 
media the LDOS approximation is excellent as long as 
Qr^n/c -C 1, which for typical Forster distances of a 
few nanometers corresponds to a frequency bandwidth Q 
of order 10w d in which the LDOS approximation can be 
made, where w d is a typical optical frequency (e.g., the 
donor emission frequency). 


For the parallel configuration near the ideal mirror, we 
find Eq. (1 (’> li . but with the integrand on the right-hand 
side replaced by 


u 2 n o 

4 Vc ^ 1 


2 

3 


- h(D x ) 


[h{2Z\) — h(Ui)\ 


(C7) 


where D\ = wir da n/c, Z\ = uj\zn/c , and U\ = 
\JD\+ 4 Zf. So we can identify both a homogeneous- 
medium and a scattering contribution between the curly 
brackets. By threefold use of the identity (IC4I) it then 
follows that ft • (Gg L ^ — Gs) • ft equals 


ft 2 n 
27r 2 w 2 c 



H(n, 


nr da. ^ 

c 


where u = \/ r da + 4 z 2 . 


c c 

(C8) 


For the perpendicular configuration near the ideal mir¬ 
ror, it can be found that the integrand of Eq. (IC3I) is 
instead replaced by the slightly longer expression 


frn 
-y—c 
47TC 



2h{2Z 1 ) + 2 


sin(2Zi) \ 
2Zi ) 


h(D i) 


-h(!h) 


4 z 2 

V r L + 4 5 - 2 


'sin(C4) 
Ui 


3h(Ui) 


,(C9) 


The frequency integral can again be performed using 
the identity (IC4l) and a standard integral of the type 
f dxa:sin(a;). The resulting expression for ft • (Gg L ^ — 
Gs) - ft, and the corresponding result (IC8[) for the paral¬ 
lel configuration are both used in Figs. [8] and [9] 
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